TOTAL NONNEGATIVITY OF INFINITE HURWITZ MATRICES 
OF ENTIRE AND MEROMORPHIC FUNCTIONS 
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Abstract. In this paper we fully describe functions generating the infinite 
totally nonnegative Hurwitz matrices. In particular, we generalize the well- 
known result by Asner and Kemperman on the total nonnegativity of the 
Hurwitz matrices of real stable polynomials. An alternative criterion for en- 
tire functions to generate a Polya frequency sequence is also obtained. The 
results are based on a connection between a factorization of totally nonnega- 
tive matrices of the Hurwitz type and the expansion of Stieltjes meromorphic 
functions into Stieltjes continued fractions (regular C-fractions with positive 
coefficients) . 
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1. Introduction 

Functions mapping the upper half-plane of the complex plane into itself (7Z- 
functions) are well studied and play a significant role in applications. The subclass S 
of 7Z- functions, the functions that are regular and nonnegative over the nonnegative 
semi-axis (also known as Stieltjes functions) is of particular interest. In this paper 
we demonstrate a connection of meromorphic 5-functions with total nonnegativity 
of corresponding Hurwitz-type matrices (Theorem 1.4). As an application, we study 
the following problem on the distribution of zeros. 

A polynomial with no roots with a positive real part is called quasi-stable. Asner 
(see [3]) established that the Hurwitz matrix of a real quasi-stable polynomial is 
totally nonnegative (although there are polynomials with totally nonnegative Hur- 
witz matrices which are not quasi-stable). A matrix is called totally nonnegative if 



2010 Mathematics Subject Classification. 30D15 (primary); 15A23, 30B70, 30C15, 34D20, 
40A15 (secondary). 

This work was financially supported by the European Research Council under the European 
Unions Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 259173. 



ALEXANDER DYACHENKO 



all of its minors are nonnegative. In addition, Kemperman (see [15]) showed that 
quasi-stable polynomials have totally nonnegative infinite Hurwitz matrices. 

It turns out that the replacement of finite Hurwitz matrices with infinite Hurwitz 
matrices allows us to prove the converse: a polynomial is quasi-stable if its infinite 
Hurwitz matrix is totally nonnegative. The key to this is given in [9]: a special 
matrix factorization, which was successfully applied to a closely related problem 
in [10]. Moreover, when a theorem involves an infinite Hurwitz matrix, it is natural 
to suggest that it can be generalized to entire functions or power series. The 
first goal of the present paper is to obtain the following extension of the results 
from [3], [15] and [9] to power series, including the converse result. 

Theorem 1.1. Given a power series f{z) = z-> X^fcLo fk zk * n ^ e complex vari- 
able z, where /o > and j is a nonnegative integer, the infinite Hurwitz matrix 



H f = 



//o fi h h h 

o h h h h 

o h h h h 

o o h h h 

o o/o/ 2 h 

\\ ; i i ; 



is totally nonnegative if and only if the series f converges to a function of the form 

/w=c ^w n "( l+ ") n -(;;-)( 1+ -) , p., 

n> (1 + i) (1 - i) 

where C, 71,72 ^ 0, x^,y\ > 0, Rea^ ^ 0, Ima,, > and 

E,t + ^Re(^) + £^+EA^<oo. 

Remark 1. Stating herein that a power series converges, by default we assume it 
to be convergent in a neighbourhood of the origin. Moreover, where it creates no 
uncertainties we use the same abbreviation for the series and a function it converges 
to. 

Remark 2. It is possible that {x^}^ n {yx}\ 7^ in the expression (1.1). If so, 
the coinciding negative zeros and poles of the function f(z) cancel each other out, 
while its positive poles remain untouched. For example, although the series J2kLo z> * 
satisfies Theorem 1.1, it converges to the function y— ^ with a unique positive pole. 
The number of such cancellations may be infinite, however it cannot affect the 
convergence of involved infinite products. 

Our second goal is achieved by Theorem 1.2, which is an extension of [10, The- 
orem 4.29]. 

Theorem 1.2. A power series f(z) — J2T=o $ kzk w ^ /o > converges to an 
entire function of the form 



/W = /oe^n( 1 



z 
a„ 



(1.2) 
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where 7 ^ 0, a v > for all v and ^2 U — < oo, if and only if the infinite matrix 
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is totally nonnegative. 

This theorem complements the following well-known criterion established by 
Aissen, Edrei, Schoenberg and Whitney. 



Theorem 1.3 ([1, 2, 6], see also [12, Section 8 §5] 
ries f{z) = ^2^ =0 fkZ k , /o > 0, the Toeplitz matrix 



T{f) 



(h 
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/o 
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.A 

/o 
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/o 



Given a formal power se- 
\ 

(1.3) 



is totally nonnegative if and only if f converges to a meromorphic function of the 
form: 

f(z)=he" ) "4, (1-4) 



n. 



1 2- 

1 Pu 



where 7^0, a Ul (5^ > for all fj,, v and ^2 V ■£- + J^ 4- < 00. 

If we require the series f(z) to represent an entire function under the assumptions 
of Theorem 1.3, we obtain that it has the form (1.2). We prove Theorems 1.1 and 1.2 
in Section 4. 

A sequence (fk)kLo is commonly called totally positive (e.g. [2]), or a Polya 
frequency sequence (e.g. [12]), whenever the matrix T(f) defined by (1.3) is totally 
nonnegative. By Theorem 1.3, the general form of its generating function is given 
by the formula (1.4). 

Definition 1. We denote by TZ (TZ~ l resp.) the class of all meromorphic 1 functions 
F(z) analytic in the complement of the real axis and such that 

Im.F(z) / ImF(z) 



Imz 



^0 



Imz 



s: for F £ Tl~ L resp 



Note that it is a straightforward consequence of the definition that 1Z- and 1Z~ 1 - 
functions are real (i.e. map the real line into itself). Furthermore, our definition 
includes real constants (like in [13]) into both classes 1Z and ~Rr x although some- 
times they are excluded in the literature (e.g. [22]). 



In general, the condition to be meromorphic is replaced by less restrictive F(z) = F(z). Basic 
properties of ^-functions can be found, for example, in [13] and (for the meromorphic case) in [22]. 
For brevity's sake we confine ourselves to meromorphic functions only. 
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Definition 2. Denote by S the subclass of ^-functions that are regular and non- 
negative over the nonnegative reals. (Since 5-functions are meromorphic, they can 
have only negative poles and nonpositive zeros.) 

Consider the infinite Hurwitz-type matrix (i.e. the matrix of the Hurwitz type) 



H(p,q) 



(bo 


hi 


b 2 
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In 


h ... 





no 


«'i 


(12 


«3 


04 
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b 2 
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64 ... 








«o 


«i 


a-2 


a 3 ... 



I : 





bo 
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b 3 ... 



(1.5) 



where p(z) — J^kLo a k zk an d q{z) = X)fc°=o bkZ k are formal power series. Given 
two arbitrary constants c and j3, we also consider the matrix 



J(c,(3) 
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13 
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/? 



V 



(1.6) 



Matrices of this type will appear in our factorizations below. 

Finally, for an infinite matrix A = (fljj)™- =1 and a fixed number p, < p ^ 1, 
we consider the matrix norm 



\\A\\ p :=supYJ/> 



.i- 1 ! 



;>i 



i=i 



Remark 3. Convergence in this norm implies entry-wise convergence. Moreover, 
the norm ||^4|| p of a matrix A coincides with the norm of the operator 



A p : 14 ^l-diag(l,p,p 2 ,. 



x, 



1,2, 



acting on the space l^, of bounded sequences. 

Remark 4. Let functions g(z), p(z), q(z) and g( k \z), p( k \z), q^ k '(z), k 
be holomorphic on D p :— {z e C : \z\ ^ p}. Then the condition 

lim \\T(gW)-T(g)\\ p =0 

K—tOO 

is equivalent to the uniform convergence of g^ >{z) to g(z) on D p , and the condition 

lim \\H(p^,q^)-H(p,q)\\ p = 
h—>oo 

is equivalent to the uniform convergence of p^ (z) to p(z) and q( k > (z) to q(z) on D p . 

Now we can formulate the more important result of this paper concerning prop- 
erties of S- functions. It has its own value apart from the proofs of Theorems 1.1 
and 1.2. 



Theorem 1.4. 



Consider the ratio F(zj — — — - 

and q(z) = X^fclo bkZ k , normalized by the equality p(0) 
conditions are equivalent: 



— - 1 — - of power series p(z) = Y^k=0 a k zk 



Of) 



1. The following 
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(i) The infinite Hurwitz-type matrix H(p,q) defined by (1.5) is totally non- 
negative. 
(ii) The matrix H(p, q) possesses the infinite factorization 

H(p,q)= lim (j(& ,/3o)J(l,/3i)---J(l,/3 J ))ff(l,l)T( 5 ) (1.7) 

converging in \\ ■ \\ p -norm for some p, < p ^ 1. Here b ^ and the 
sequence (/3j)j>o * s nonnegative, has a finite sum and contains no zeros 
followed by a nonzero entry, that is 

POjPi,. • • ,Pw-i > 0, p w = p w +i = • • • = 0, 

°° (1 8) 

^ d ^ oo, and \, Pi < °°- 

3=0 

The matrix T(g) denotes a totally nonnegative Toeplitz matrix of the 
form (1.3) with ones on its main diagonal. 
(iii) The ratio F(z) is a meromorphic S -function; its numerator q{z) and de- 
nominator p(z) are entire functions of genus up to a common meromor- 
phic factor g(z) of the form (1-4), g(0) = 1. 

Remark 5. Note that 

J(c,0)H{l,l)=H(l,c). (1.9) 

If w is a finite number in (1.8), then /3 u -\-i = fi^+2 = ■ ■ • = implying 

J(b , p" ) • • • .7(1, Pu+i) 77(1, 1) - J(b , po) ■ ■ ■ J(l, /3 w+ i) JfrPu+a) 77(1, 1) = • • • . 
As a consequence, the factorization (1.7) can be expressed as follows in this case 

H( ,ij(bo,0)H(l,l)T(g) = H(l,b o )ng) if« = 0; 

[P,Q) \j(b 0l l3 )J(l,/3 1 )---J(l,^_ 1 )H(l,l)T(g) if0<w<oo. 

Remark 6. The number p in Theorem 1.4 can be anywhere in (0, 1] n (0, po), here 
po denotes the radius of convergence of g(z) (which is positive by Theorem 1.3). 

The matrix T(g) from the condition (ii) of Theorem 1.4 is the Toeplitz matrix 
of the function g(z) from (iii) given by (1.5). 

If we require p(z) and q(z) to be entire functions in Theorem 1.4, then the 
function g(z) has the form (1.2) or g{z) = 1 and (1.7) converges in || ■ ||i. 

Remark 7. In the case g(0) = 6o = it can be convenient to "trim" the ma- 
trix H(p, q) by removing its first row and its trivial first column. This corre- 
sponds to replacing J(0,po) in the factorization (3.9) by its diagonal analogue 
diag(l,pb,l,pb,..-)- 

Remark 8. Since entire functions of genus have unique Weierstrafi' representa- 
tion, it makes sense to consider the greater common divisor of a subset of this 
class. Accordingly, two entire functions p and q of genus are coprime whenever 
gcd(p,q) = 1. 

Consider the continued fraction 

b + -^- + -±- + -^ + + ^_i£, 6 0>0 , po,Pi,...,£ w -i>0, 0<«<oo, 

(1-11) 
where we combine both finite (terminating) and infinite cases. If the continued 
fraction is infinite, we assume uj = oo. The following Corollary (see its proof in 
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Subsection 2.2) allows us to connect the factorization (1.7) with continued fractions 
of this type. 

Corollary 1.5. Let F(z) = ^pfr be a meromorphic S-function, where the entire 
functions p(z) and q(z) are of genus 0. Then it can be expanded into a uniformly 
convergent continued fraction of the form (1.11) with exactly the same coefficients 
bo and (j3j)^Z , X^=o Pi < °°> as ^ n ^ le factorization (3.9) of the matrix H(p, q). 
No other continued fractions of the form 



F(z) = co + 



c 2 z r2 C 3 Z r3 c^z 7 



1 + 1 + 1 +■■■+ 1 

where Cj ^ and rj 6 N for j = 1, . . . , u, ^ co ^ oo, can correspond to the Taylor 
series of F{z). 

Remark 9. Corollary 1.5 implies that each pair {p{z),q{z)) satisfying Theorem 1.4 
determine a unique factorization of the form (1.7). 

Let p(z) and q(z) be real polynomials. Denote 

u(z) := J2 akZ n ~ k = z n p f - J and v{z) := ^ b k z n - k = z n q f- 

where n — max{degp, degg}. In this case it is more common to work with the 
matrix H (u, v) := H (p, q) instead of H (u, v). 

In fact, Theorem 1.4 extends the following result by Holtz and Tyaglov to mero- 
morphic functions. In [10, Theorems 1.46 and 3.43, Corollaries 3.41-3.42] they 
established that the matrix H (u, v) is totally nonnegative if and only if it can be 
factored as follows 

H(u,v) = J(c ,l)...J(c j ,l)H(l,0)T(g), Cl ,..., Cj >0, (1.12) 

where T(g) is totally nonnegative and g = gcd(u,w). Note that the factoriza- 
tion (1.12) corresponds to (1.10) after the substitutions 60 = c o, Po = (ci)~ 1 and 
f3i-i = (ci-iCi)^ 1 for i — 2, . . . , j. Moreover, by Theorem 3.44 from [10] the ma- 
trix H(u,v) is totally nonnegative if and only if v(z) and u(z) have no positive 
zeros and ^ e "Rr 1 . Since 

ff(*) «(i)' 

we obtain the polynomial analogue of Theorem 1.4. 

Earlier, Holtz (see [9]) found that the infinite Hurwitz matrix of a stable poly- 
nomial (i.e. a polynomial with no roots with nonnegative real part) has the fac- 
torization (1.12) with T(g) equal to the identity matrix. Additionally, each of the 
factors J(cj,l) corresponds to a step of the Routh scheme. These factorizations 
coincide because the problems considered in [9] and [10] are closely connected (see, 
for example, the monographs of Gantmakher [7, Ch. XV] and Wall [21, Chapters 
IX and X]). In order to deduce Theorem 1.1 from Theorem 1.4, we are using the 
same underlying connection. 

2. Basic facts 

Here we consider some facts that are quite significant, although, in fact, they 
are not new. We put them here to introduce the area and our notation. The 
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most "non-standard" assertion here is Lemma 2.11, since it reverses the approach 
of Theorem 2.10. 

2.1. 5-functions in terms of Hurwitz-type matrices. Consider power series 

oo oo 

p(z)=^a k Z k , o = l, and q{z)=^b k z k , o Js 0. (2.1) 

fc=0 fc=0 

Let us introduce the following notations 
p Q (z):=p(z), p-i(z) := q(z), H := H(p,q), and H := if(p ,P-i) = H. 

Denote the minor of a matrix A with rows i\,i2, ■ ■ ■ ,ik and columns ji,J2, ■ ■ ■ ,jk 

by 

Jn i-2 ■■■ 1 1. 



K 3i 32 ■■■ 3k 
In addition set 

a«-a(» J ::: *). (2.2) 

If the number /3 = 6i — 6o «i = -Hq ^ s nonzer o, we define 

g(z) -6 p(z) 

PH Z ) := 5 and ff i ;= H (PiiPo)- 

p z 

Now we can perform the same manipulations with the pair pi(z), po(z). That is, 
we can make the next step of the following algorithm. 

At the jth step, j = 0, 1, 2 . . . , the series Pj(z) and pj_i(z) are already defined, 
as well as the matrix Hj = H(pj,pj-i). We set 



/3 j:=J fff, (2.3) 



and, if j3j is nonzero, we set 



p J+1 (z) := P] -^ Z) PplMlM (note that Pj ^(0) = 1 when j > l) (2.4) 

so that Hj + i := H(pj + i,pj). These steps can be repeated unless j3j — 0. In 
Corollary 2.7 we will show that j3j > whenever Fj(z) = V3l / z V represents a 
non-constant meromorphic 5-function. To do this we need some auxiliary facts. 

Suppose that ft ^ 0, i = 0,1,..., j for some nonnegative j, such that the 
power series pj-\(z), Pj(z) and pj+\{z) arc defined according to the recurrence 
formula (2.4). 

Lemma 2.1. The identity^ 

(2 3 ... k k + l\_M\ (2 3 ... fc-1 k\ 

Hj \2 3 ... fc i + l)- P i Hj+1 \2 3 ... fc - 1 »J ' 

ZioZds /or all k ~ 2,3, . . . and i — k, k + 1, . . . . 



The notation [aj stands for the maximal integer not exceeding a. 
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Proof. Without loss of generality we consider the case j = 0, since for higher values 
of j the relations (2.3)-(2.4) are analogous. In the case k = 2m 



Po n Hi 



2 3 
2 3 



2m — 1 2m 
2m — 1 i 



= /3 m #i 



1 2 
1 2 



2m 
2m 



a a\ a 2 

fei — 6 ai »2 - 6002 
ao ai 

b\ — 6qOi 





a a,\ 
b b\ 



O-lm-1 
&2m-2 ^ 0002771—2 

G2m-3 
^ra-3 — boa,2m-3 



- 1 2m 

-1 i 

a-j-i 

Oj-l — boCLi-i 

aj-2 

Oi-2 ^ bodi-2 



« 2 
02 



ao ai 
&o 61 







02m-2 
&2m-2 

fl2m-3 
02m-3 



y-m-i - 

Oi_l 
bi-1 

a.1-2 

bi-2 



boa m -i 



Ho 



bi-r 



2 3 
2 3 



60 ai 



2m 2m + 1 
2m i + l 



here the equality ao = 1 has been used. For k = 2m — 1 the transformation remains 



the same. 

In particular, if we suppose that /3o,/3i, 
implies 



□ 



. ,/3fe_i > for k ^ 3, this lemma 



H 



(fc) 



^J F (fe-D_ 



'i+i 






= ... =H 



( 3 ) nVVJ 

fc-3 LLPi+j-1 



fc-2 



LVJ 



(2.5) 



The next theorem was established by Chebotarev, see [4] and [5, Ch.V §1]; see 
also the proof of M. Schiffcr and V. Bargmann in [22, II. 8]. At the same time, it 
can be derived as a particular case from Nevanlinna's theory, see [14, Theorem 8]. 

Theorem 2.2 ([4, 5, 22, 14]). A real meromorphic function F(z) regular at the 
origin is an IZ-function if and only if it has the form 

F(z) = B + B lZ + V {-^--^L 

^ (2.6) 



o-l 



where 2, — \~~ "^ °°' ^i > and A v < 0, a v <E K /or 1/ = 1, 2 . . . , w. 



The proof of this theorem relies on the following fact which we will use later. 

Lemma 2.3 (see e.g. [5, Ch.VI §8]). Let entire functions q(z) and p(z) have no 
common zeros and such that F = - € 1Z is not a constant. Then F' ' (z) > on the 
real line, the zeros of p(z) and q(z) are real, simple and interlacing. 

The interlacing property means that between each two consequent zeros of p{z) 
there exists a unique root of q(z) and vice versa. The proof from [5] is based on the 
behaviour of meromorphic 1Z- functions in neighbourhoods of its zeros and poles. For 
completeness, we deduce this lemma here from the partial fraction expansion (2.6). 
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Proof. Let F(z) = ^A have the form (2.6). If z is not real, then 

Imz j—t \z + a v \ z 

Therefore, F(z) (as well as q(z)) has no zeros outside the real axis. 

Now from (2.6) it follows that F{z) is real and can only have simple poles. Since 

F'(z) = B 1+ V ~ A " >0, zel, 

the function -F(z) grows between any of its two subsequent poles z\ and Z2 from 
— oo to +oo. So there is one and only one z* € (j?i, 22) such that F(z*) = q(z*) = 0. 
For the same reason, there exists a unique zero of p(z) between any two subsequent 
zeros of q(z). □ 

The next theorem is a consequence of Grommer's theorem (see [8, §14, Satz III]) 
and Theorem 2.2. It can be proved by applying the Hurwitz transformation [11] 
(see also [14, §6.1], [5, Ch.I §7], [10, Theorem 1.5], [7]) to the matrices of the Hankel 
forms corresponding to F(z). 

Theorem 2.4 (e.g. [5, Ch.V §3]). A meromorphic function F(z) = ^f\, where 
p(z) and q(z) are of the form (2.1), is an IZ-function if and only if there exists I, 
^ I ^ oo, such that 

H (2m+i) > 0, m = 1, 2, . . . , /, H( 2l +V = H^+V = • ■ ■ = 0. 

Moreover, I is finite if and only if F(z) is a rational function with exactly I poles, 
counting a pole at infinity (if exists). 

Let f3i y^ 0, i — 0,1, . . . ,j for some nonnegative j, and the power series pj~i(z), 
Pj(z) and pj+i(z) be defined by the recurrence formula (2.4). Suppose that the 
ratio Fj(z) = P3 ~) V of formal power series converges to a meromorphic function. 
Then there exist entire functions pj-i(z) and pj(z) with no common zeros such 
that 

F J ^) = E ^T$, Pi-i(0)=W-i(0) and p,(0) = Pj (0) = 1. 
Define the power series g(z) := v ~ fl satisfying g(0) = 1. Then 

- I \ Pj( Z ) J - / \ Pj-l( Z ) 

If Fj € 7?. is a non-constant function, then by Theorem 2.4 the inequality j3j — 
H-'>0 is satisfied. So from the formula (2.4) we find 

*j+l{Z) - 



Pj+i( z ) Fj(z) ~ P3-1W 
Lemma 2.5. The ratio Pj+1 z converges to the entire function 

^(^ftdWlMftW. (2.7) 

TTie pairs (pj(z), j»j+i(z)) and (p J _ 1 (z),pj +1 (z)) have no common zeros. 
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Proof. Dividing (2.4) by g(z) gives ^/A — p*j+i(z), that means the relation (2.7) 
holds. Consequently, 

Pj-i(z) = 0jZPj+i(z)+Pj-i(O)Pj(z)- 

Each common zero of any two summands in this equation must be a zero of the 
third summand. Since the functions pj-i(z) and pj(z) have no common zeros, the 
pairs (pj^i(z) 1 pj + i(z)) and (pj{z),Pj+i(z)) also have no common zeros. □ 

This lemma implies that Fj+i(z) represents the meromorphic function 

Pj( z ) 



F j+ i(z) 



P] + l( z )' 



Lemma 2.6. If the meromorphic function FAz) is not a constant, then Fj £ S if 
and only if Fj,Fj + i £ TZ and Fj(0) ^ 0. 

Proof. Let Fj £ S, then Theorem 2.2 gives that it has the form 
F j (z) = B +B 1 z+ £ f^--^-)=B Q + B 1 z + z £ t^M , 

where V" —^ < oo, B ^ 0, Bi > and A„ < 0, a u > for v = 1, 2 . . . , u. 

It is enough to show that Fj+i(z) is a well-defined 7vL-function. Consider the func- 
tion 

Gj(z):=^tll^m =Bl+ V (-^K) = Bl + V (AlM. 

— z ^-^ —z + a v ^-^ z — a v 

It has the form (2.6) and, hence, is a meromorphic 7?.-function by Theorem 2.2. 

The mappings z 4 ^ and z h-> ~z are in the class TZ~ l {i.e. they map the 
upper half of the complex plane into the lower half of the complex plane). Since 

/3j — Hj > and Gj £ TZ, the function composition 

- ° G > ° <- ' )) W = rjS = F 1 \ — (M = F >+^ 
J Gj{-z) Fj(z) - pj-i(0) 

is an 7?.-function as well. 

Conversely, let Fj,Fj+i £ TZ and Fj(0) ^ 0. The inequality f3j > holds, 
therefore i^+i (z) ^ and the meromorphic function 



ft F 3 -(-z)-F 3 -(0) 

'""- F j+1 {-z)- 



Gj[z) ■-- 



is an TZ- function. On one hand, Theorem 2.2 gives 

F j (z)=B + B l z + z J2 { ~ Av,<7v \ where 



Z + Vv 



l<I/<OJ 



V ^—4- < oo, Bi > and A v < 0, av € R for i/ = 1, 2 . 
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such that 

F J (-z)-FM _ _ ^ {A v /a v ) 



G s {z)= r "-<>_-/> w =B 1+ £ 



z — a v 

1<i/<o; 

On the other hand, Theorem 2.2 states that each 7\L-function has negative residues 
at its poles. That is, ^- < for all v since Gj G 7Z. Therefore, the poles — er„, 
v = 1,2. . . ,u>, of the function Fj are negative. Consequently, Fj £ S. D 

Corollary 2.7. Suppose that for some j ^ t/ie function Fj(z) is in the class S. 
Then ft Js 0. XTie inequality ft > implies F-j+i G 5, w/izfe i/ie equality ft = 
implies that Fj(z) is constant. 

Proof. If Fj(z) is a constant then ft — H } = by Theorem 2.4. Let Fj(z) be a 
non-constant 5-function. Applying Lemma 2.6 to it gives Fj+i G 7?.. Consequently, 
ft+i = if Fj+\{z) is a constant and ft+i > if it is not. Moreover, we have 
Fj+i (0) = 1 > 0. Thus, the corollary holds in the cases of constant Fj (z) or Fj+i (z). 
Suppose that ft,ft + i > 0. Then Lemma 2.1 implies that 

H (2 m+ 3) = p\^\ [^\ H {2rn + l) = pm+1 prn_ ^f™^ , m = 1, 2, . . . . (2.8) 

That is, for each natural m the sign of H _- + ™ coincides with the sign of iJ- m . 
Since F,- G ft, Theorem 2.4 yields F i+2 G K. That is, F i+1 G S by Lemma 2.6. □ 

Theorem 2.8. ^4 meromorphic function F(z) = ^f), where p(z) and q(z) are 
series of the form (2.1), is an S -function if and only if there exists uj, 2 ^J iq ^ oo, 
such that 

H (k) >0, fc = 2,3,...,w and H^ +1) = H^ +2) = ■ ■ ■ = 0. (2.9) 

Moreover, u> is finite if and only if F(z) is a rational function with exactly I ^^- J 
poles, counting a pole at infinity (if exists). 

Proof. By definition, Hq = H^ = 1 > 0. Denote po(z) '■= p(z) andp_i(z) := q(z) 
such that F(z) = Fq(z). From the recurrence formulae (2.3)-(2.4) we obtain the se- 
quences (Pj)j=-i an( i (ft)^=o 2 > wnere ft ^ for all j = 0, 1, ... w— 3 and 2 ^ w ^ oo. 
Whenever cj < oo we also have ft>-2 = 0. 

Suppose that F G S. Then ft > for all j = 0, 1, . . .w - 3 by Corollary 2.7. 
Furthermore, the identity (2.5) gives 

ff^-ff^-n^-i >0 ' i = 3,4,...w. (2.10) 

Let w < oo, then F u —2(z) is a constant by Corollary 2.7. Therefore, we have 
H^_ 2 = H^_ 2 — ■ ■ ■ = since all these minors contain proportional rows. By the 
identity (2.5), this is equivalent to H^ u+r > = H iuJ+2 *> = • • • = 0. 

So we obtained that F G S implies (2.9). The number of poles the function F(z) 
has can be determined from Theorem 2.4. 

Now suppose that the conditions (2.9) hold. If w = 2 then H^ = H^ = ■ ■ • = 
and the assertion of this theorem is equivalent to Theorem 2.4. In the case of 



-j - 1, and 
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3 sj uj ^ oo we have /3o > 0, so by Lemma 2.1, 

^(2m+l) =/3 -L^ 1 J jff (2m+2) > ^ m = 1, 2 , . . 

Hence, the functions -F(z) and Fi(z) are 7?.- functions by Theorem 2.4, and 
Lemma 2.6 yields F £ S. □ 

2.2. 5-functions as continued fractions. A continued fraction of the form 

Ciz'" 1 C 2 Z r2 C 3 Z r3 C w 2 r ^ 

F( 2) = Co + n _ + n _ + n _^ + ^_, where (jm) 

Cj 7^ and r 3 eN for j = 1, . . . ,u, ^ u> ^ oo, 

is called a (general) C -fraction. The special case of (2.11) that corresponds to 
rj = 1 for all j = 1, ... ,u> is called a regular C -fraction. Continued fractions of the 
form (2.11) are able to represent power series uniquely, that is the following fact is 
true. 

Theorem 2.9 ([16], see also [18, §21, Satze 3.2-3.5, 3.24]). Each (formal) power 
series F(z) — J2T=o Skzk corresponds to a fraction of the form (2.11). This corre- 
spondence is set by the following sequence of relations 

F (z) = F(z), Co=F(0), Fi(z)= ^ i = 1,2, ... ,uj, (2.12) 

where uj ^ oo is such that Fi-i(z) ^ _Fi_i(0) for i— 1 < uj and F ul {z) = F u (0). The 
exponents ri are positive integers chosen together with the complex constants Ci in 
such a way that Fi(0) = 1. 

// two C -fractions (finite or infinite) of the form (2.11) correspond to the same 
power series, then they coincide. A C -fraction is finite if and only if it corresponds 
to a rational function (and, hence, represents that function). 

Moreover, if an infinite continued fraction of the form (2.11) converges uni- 
formly in a closed region T containing the origin in its interior, it represents a 
regular analytic non-rational function of z throughout the interior of T. Further, 
the corresponding power series converges to the same function in and on the bound- 
ary of the largest circle which can be drawn with its center at the origin, lying wholly 
within T . 

Suppose that F(z) = ^A is an S- function. We again denote Fq{z) :— F(z), 
Po(z) :— p{z) and p^\(z) :— q(z) and use the recurrence formulae (2.3)-(2.4) to 
obtain the sequences (pj) ( j = _ 1 and (/?;,- )^ =0 , where f3j ^ for all j = 0, 1, . . .ui — 1 
and — 1 ^ uj ^ oo. In the case uj < oo we also have j3 u = 0. 

For each j = 0, 1, . . .u> — 1, we apply Corollary 2.7, obtaining j3j > and Fj € S. 
If uj is a finite number, then F u is a constant. From the relation (2.4) we have 

P ^ = ^ = ft-^ + A' i = 0,l,..., W -l. (2.13) 

Pj{Z) %l(2J 

These formulae can be combined into the continued fractions 

F(z) = F (z) = b + — + — + — ++ —j— and (2.14) 



PjZ Pj+lZ f3j+2Z Pu-l Z 

1+1 + 1 +•••+ 1 



Fj{z) = 1 + lll 2±±1 t^l 20=11, (2.15) 
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where j = 0,1, ... ,lu — 1 and /3i > for i = 0, 1, . . . , u — 1. These are regular 
C-fractions, and the relations (2.13) set the correspondence satisfying (2.12). That 
is, the continued fractions in (2.14) and (2.15) corresponds to F(z) and Fj(z) for 
all j , respectively, by Theorem 2.9. In particular, they are finite if and only if F(z) 
is rational. 

Furthermore, there is a power series g(z) such that pj(z) := P3 ,V are entire 
functions for j = — 1, 0, . . . , uj, and for j ^ pj~i(z) and pj(z) have no common 
zeros (see Lemma 2.5). Observe that the relations (2.13), (2.14) and (2.15) remain 
the same, if we replace all the series Pj(z) by the functions Pj{z). 

It is convenient to study the continued fractions (2.14) and (2.15), using the 
following. 

Theorem 2.10 (Stieltjes, [20, n os 68-69]; see also [18, 21] 3 ). Let b Q ^ 0. A sequence 
of positive numbers /3q, /3i, . . . , /3 W , — I^uj^oo, has a finite sum if and only if the 
continued fraction (2.14) converges to a meromorphic S -function and its partial 
numerators and denominators converge to coprime 4 entire functions of genus 0. 
That is, if the jth convergent (approximant) to F(z) is denoted by p. (I , then for 
j — > oo we have 

Pj(z) — » p(z), Qj(z) — > q(z) and 
(3 z (3 lZ fij-iz Qj{z) q(z) 

k) + — + — + - + ^ = T^z) ^ pjz)= F{Z) ' 

where p(z) and q{z) are coprime entire functions of genus 0. The convergence is 
uniform on compact subsets of C containing no poles of the function F{z). 

To apply this theorem we need to distinguish the case of J^ ■ /3j < oo 

for uj = oo. 

Lemma 2.11. Let the functions p\(z) and po(z) be entire of genus 0, coprime and 
such that their ratio S 3 F\ := S 2 - is not rational. Then pj(z) — > 1 as j — > oo 
uniformly on compact subsets ofC, and J27Li Pj < oo. 

Proof. According to (2.7), all the functions Pj(z), j = 0, 1, . . . , are entire of genus 
and pj(0) = 1. Moreover, po, ■ ■ ■ ,pj are coprime (by Lemma 2.5) and hence have 



representation is valid for j = 0,1 



only negative zeros (since by Corollary 2.7 Fj = ^L-i. g £). Therefore, the following 



CO OO • \ 

^ ) =^ a ov=n( i +^))' 

where < a[ ^ o^ ^ • ■ • and Y^,vLi ~UT < °°- The coefficients ai , k = 1, 2, . . . , 

er, y 

are equal to 

oo oo oo 1 



U ' ' ' 



= 1 »a=l ik = l CT n a i 2 " ' a ik 



J The separate convergence of numerators and denominators was shown by Slcszyriski in [19] . 
This fact was obtained by Maillet in [17]; see also [18, p. 150]. 
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Note that these sums are convergent, since 5 

oo oo 



a? 



<EE-"E— — ■•• — = («i j) ) whenfc = 2,3,4,.... (2.17) 



i 1= li 2 = l J fc = l U i x "j 2 °i,, 



By Lemma 2.3 the zeros of pj (z) and pj-i (z) (which are negative) must be simple 
and interlacing. In addition, F'Az) > for real z implying that < of < erf . 
Hence 

< a[ j ' 1] < af ] < a^ < o£> < o£ _1) < v® < . . . . (2 .i 8 ) 

Now we estimate the expression (2.16) using the inequalities (2.18) and obtain 

( 7 ) ( 7 — 1 ) (1) 

that ^ af < af for j = 1,2, . . . , i.e. the sequence of positive numbers af 

(i) 
decreases in j for fixed k. Therefore, there exists a finite lim J _s. 0O af ^ dependent 

on k. 

At the same time, the equality (2.7) implies that the first Taylor coefficient af 

for any j has the form 



Consequently, 



V ft + af' and V ft = a?' - lim af'. (2.19) 

i=l j=l 



Therefore, the series Yl'jLo Pj converges. 

As a consequence, for an arbitrary positive number R there exists an integer 
jo(R) such that fti? < 4 for all j ^ jo- By virtue of Worpitzky's test (as it stated 
in [21, p. 45], see also [23]) the continued fraction 



F j0 (z) = l + 



PjpZ ft +l^ /3j„+2 z 

1+1 + 1 



converges to an analytic function uniformly in the disk \z\ < R. This analytic 
function coincides with Fj (z) (since Fj (z) corresponds to the continued fraction, 
see (2.15)). Therefore, pj has no zeros in this disk, that is R < of"' < of , j ^ jo- 
Letting R tend to infinity, we obtain lim^oo of = oo. According to (2.16) we 
have 

oo 

a\" = 



l l ~ Z^ (J) 



-1 O, 



s=l 



5 In fact, even an estimate stronger than (2.17) is valid (cf. [19, p. 105]). For each tuple of 
itinct numbers (ii , «2 , • • • , ifc) there is onl; 
side of (2.16). At the same time, the sum 



distinct numbers (ii, 12, ■■■ ,ik) there is only one summand [erf <rf ■ ■ ■ erf ) in the right-hand 



OO OO OO _. -. ,. 

il =1*2 = 1 »fc = l ff ii a io ^i,. 



contains exactly A;! summands of this form. Therefore, 

1 

ki 



af ] < — (af^ ) for k = 2, 3, 4, 
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and each term in this series monotonically tends to zero as j — > oo. For any e > 
there exists N such that 

OO _. OO -. oo ^ 

E — < £ ' which for 3 > 3o gives J2 —^ E — <£ - 



i=N+l u i i=N+l u i i=N+l u i 

N 



On the other hand, lim > —7^- = 0, so the coefficient a\ also vanishes. Now 

from (2.17) we obtain that pj(z) —¥ 1 as j — > 00 uniformly on compact subsets 

of C. D 

Corollary 2.12. Under the assumptions of Lemma 2.11 there exists a positive 
number M independent of j such that 

WHfa+uPJWi < M for j = 0,1,..., 

where the matrix H(pj+i,pj) is defined by (1.5), and 

\\H(p 3+1 ,p 3 ) - #(1,1)11x^0. 

Proof. Since ag ^ ajr ^ for all j, fc = 0, 1, ... , we have 

{OO OO ■ OO OO 

E«* +1) ,E% a) =E4' ) <E4 0) =po(i)<oo. 
fe=0 fc=0 J fc=0 fc=0 

Observe that aj) = 1 whenever j ^ 0. Therefore, by Lemma 2.11 we obtain the 
required 

l|i?(Pi+i,Pi) - #(M)lli = max<j E°* +1) 'E°* 



0+1) y^o) 
fc=i fc=i 



OO 

2a«-ft(l)-li^0. 

fe=i 

a 

Corollary 1.5. Let -F(0) = ^r\ be a meromorphic S -function, where the entire 
functions p(z) and q(z) are of genus 0. Then it can be expanded into a uniformly 
convergent continued fraction of the form (2.14) with the coefficients bo and (f3j)^Z , 
12^=0 Pj < °°; given by (2.3). No other continued fractions of the form (2.11) can 
correspond to the Taylor series of F(z). 

Proof. If F = 2 g S then F(z) can be formally developed into the contin- 
ued fraction (2.14) with the coefficients bo — F{0) and (/3j)£T given by (2.3). 
By Lemma 2.11, the coefficients of this continued fraction satisfy the condi- 
tion Y^j=o @j < oo. Consequently, Theorem 2.10 implies that (2.14) converges 
uniformly on compact sets containing no poles of its limiting function. So the rest 
of the proof comes from Theorem 2.9: the continued fraction (2.14) corresponds to 
and converges to F(z), and there is no other continued fraction of the form (2.11) 
corresponding to F(z). □ 
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3. Total nonnegativity of Hurwitz-type matrices 

This section contains the proof of Theorem 1.4 and Corollary 1.5, preceded by 
several auxiliary facts. 

Suppose that meromorphic functions p and q are regular at the origin and have 
the Taylor expansion (2.1). Consider the Hurwitz-type matrix H(p,q), defined 
by (1.5). 

Lemma 3.1. If - G S and there exists a meromorphic function g(z), such that 
the ratios p(z) := £pi and q{z) := $¥l are entire of genus and coprime, then the 

matrix H(p, q) can be factored as in (1.7), where the numbers Pj, j = 0, 1, . . . , are 
given by (2.3) (possibly followed by zeros). Moreover, 

\\H(p, q) - J(b ,j3 ) J(l, fa) ■ ■ ■ J(l, Pj) H(l, 1) T(g)\\ p ^^ 0, 

where p, < p ^ 1, is such that g(z) has no poles in the disk \z\ ^ p. 

Proof. For any two matrices A = (afe/)^°/ =1 and B = (&fc/)£°; =1 such that ||^4||i < oo 
and || -B ||p < oo the following estimate (implying the existence of the product AB) 
is true 



oo > |L4||i||£|| 



sup yjctoi sup y^i&n 

lS^£;<oc , 1 l?Cm<oo . , 



y'- 1 > 



sup ^lafezl^l^V X > SU P Yl 



2 , a klblj 



1=1 



nJ- 1 



\AB\ 



(3.1) 



(3.2) 



=i j=i 

Now we note that the decomposition 

H(p,q) = H(p,q)T(g) 

is valid. It can be checked by the straightforward multiplication. 

Denote po(z) := p(z) and p~i(z) ■— q(z). We are now using the algorithm (2.3)- 
(2.4) to construct the (longest possible) sequence (pj)^__ 1 of entire functions, 
^ u) ^ oo. By Corollary 2.7, the corresponding numbers (pj)^_ satisfy Pj > for 
all j = 0, 1, . . .u> — 1. In the case of finite uj we have Puj-i(z) = p w -i(0), Puj(z) = 1 
and p u = 0; we extend the latter equality by P^+i = Pu+2 = ■ ■ • = 0. 

The identity (3.2) implies the factorization (1.10) in the case of q(z) = q(0)p(z) 
(corresponding to w = 0). Suppose that uj > 0. If we expand Pi(z) and pi_i(z), 
i = 0, 1, . . .j < uj, as follows 



P 



(z) = ^2c k z k and Pi-i{z) = ^ d k z k , 



k=0 



fc=0 



then the matrix H 



t+i •= 



H{pi + i,pi) takes the form 



H, 



i+i 



/co Ci c 2 c 3 

-j-{di-d ci) j-{d 2 -di ) C2) ^-(cod 3 - d c 3 ) 

c ci c 2 

■fiidi - d ci) ■j-{d 2 -d C2) 

Co ci 
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Left-multiplying this matrix by 



./(do, A) 



(<k 


ft 











... 








1 








... 








do 


ft 





... 














1 


... 




\ : 











do 


ft ... 



we obtain Hi. On putting i successively equal to 0, 1, 
find that 



. , j and applying (3.2) we 



H(p i q) = J(b o ,0o)J(l,0 1 )---J(l,0 j )H(p j+lt p j )T(g). 



(3-3) 



The finite product of the matrices is well defined and associative, since the matrices 
■7(1, •) and H(l, 1) have at most two nonzero entries in each row and column. So 
if uj is a finite number, then for j = u; — 1 the equality (3.3) coincides with (1.10). 
Since ||T(g)|| p = g(p) < oo, from (1.9), (3.1) and Corollary 2.12 we obtain the 
assertion of the theorem for uj < oo. 

Suppose that lu = oo and let us prove that the difference between the product 
in (3.3) and the right-hand side of (1.7) converges to zero as j — > oo. There exists 
an index jo ^ 1 such that 



5>4 



(3.4) 



Let 

V:= J(6o,ft)J(l,ft)---J(l,ft -i) and Uj := J(l,ft ) • • • J(l,ft) 
for j = jo, jo + 1, • • • ■ Then we can express the equality (3.3) as follows 

H(p,q)=VU j H(p j+1 ,p j )T(g). 

The matrix Uj = J(l,ft, ) • • • J(l,ft) is upper triangular and has no negative 
entries since it is a product of upper triangular matrices with nonnegative entries. 
The diagonal elements of Uj are the products of corresponding diagonal elements 
of J(l,(3j ), . . . , J(l,/3j) and, thus, are equal to 1 on the odd rows and on the 
even ones. 

More specifically, denote the entries of Uj by v/fJ so that 



U, 



(i) 



fej=i 



For j Js jo, k,m = 1,2,... the equality Uj+\ = Uj J(l, ft+i) implies 



a k,l ~ "fe,l' 



a k.2m ~ a k,2m~l h J J+l> 



The following entries for all j ^ jo must be zero 



.0) 



,(i) 



0. 



1,2, 



a k,2m+l ~ u k,2r, 



2m, 2m + 1, 



l k,2r. 



u, k,2m-\ ~ u 'k,2m 

For j = jo we have Uj = J(l, ft), so the nonzero entries of [/, are only 



,0') 



,« 



,W 



= fti u 2m,2m+l = 1 and «2m+l,2m+l = !> where m=l,2, 
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Consequently, the following estimate is valid 

(i \ m -^ 
«Mm + «£L+i< (Eft ' m=1 » 2 .---. fc = l,2,...,2m. (3.5) 

Suppose that (3.5) holds for some j ^ jo, then for fc ^ 2m we have 
,0'+i) , „,0'+i) _„,0') 



-i ft-+i + (4!L + «S m +i) < 



$: 



j \ / j X^-LU / i+1 \— L4J 

ft+i + E^ Eft < Eft 

i= Jo J \i=jo ) \i=Jo 

By induction, the conditions (3.5) hold for all j ^ jo- Therefore, by (3.4), 

U U) +U U) < /y-oo fl \"»-L4J <o-" l +L*J 

where m = 1, 2, . . . and k = 1, 2, . . . , 2m. As a consequence, 



/ oo \ oo 

1^11!= sup ug + EK'L+4 j ;L + i)HE 2 " 



-m 

m— 1 / m— 



Since uj^ - ug = and u^_i - «£L_i = «*,L-2> "i > 1, we have 



oo 

E/1 (i+i) 

l&fc<oo , vl 



\u i+i -u s h = . sup 2. 1K,^-i -^L-J + R'S -«Sm 

'm=l 



< 



/ OO oo oo \ 

-p E »2L- 9 + E «2L + E «ffi? < 

i<fc<oo \ m=2 m=1 m=1 y 

(oo oo \ 

2 E «2L + E «g£ < 
m=l m=l / 

(OO oo \ 

2ft" E «&m-l + ft+1 E ^L-l < 
^1 m^l / 

2ft- 1| J7,_i ||i + &+i||17,||i < 4ft- + 2& +1 ^^> 0. 



That is, (Uj)°°=j is a Cauchy sequence, hence it converges to its entry- wise limit [/*. 
So we have shown that ||t/j||i is bounded uniformly in j and there exists U* satis- 

fymg||^-LT*||i^^0. 

Using (3.1) we obtain ||V||i < || J(&o,A>)||i|| J(l,j8i)||i • • • II J(l,& -i)||i < oo 
and 

\\VU j H j+l T(g)-VU,H(l,l)T(g)\\ p ^ 

\\V\\i \\Uj H j+1 - U, H j+1 + U* H j+1 - U* H(l, 1)^ \\T(g)\\ p < 
II^IIJI^-^llill^+ilUlr^ll. + II^IIJ^^JI^+x-^l.^llJlT^II^ 

(3.6) 
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The expansion of g{z) into a power series at the origin converges for \z\ ^ p abso- 
lutely, hence ||T(g)|| p < oo. According to Corollary 2.12, 

||J? i+ i||i = \\H{pj + x,Pi)\\x < M and \\H(p J+uPj ) - #(l,l)||i ^^> 0, 
so the right-hand side of (3.6) vanishes and, therefore, 
H(p,q) = VU*H(l,l)T(g), where VU* = Urn (J(6 ,A)) .7(1, ft) •• • J(l,&)) ■ 



Remark 10. Applying this lemma to the ratio 
the matrix [/*. Since 



n 



J0 7 , we can explicitly determine 

PjQ \ Z ) 



,(i+i) 

l k,2m 



/3,+i^L-i^^O form =1,2,. 



and k = 2m, 1m + 1, 



from Hj = U*H (1,1) we get 



^ = F ioJ ff T (0,l), 

where A T stands for the transpose of a matrix A. 

Now consider mcromorphic functions p(z) =: po( z ) an d g(z) =: p-\(z) ^ 0, 
(/(O) ^ 0, with the power series expansions given by (2.1). Suppose that /3o, /3i, 
. . . , /3j-i 7^ 0. Then we can define the functions pi(z), P2(z), ■ ■ ■ , Pj{ z ) via the 
formula (2.4). 

Lemma 3.2. If the Hurwitz-type matrix Hj := H(j)j,pj^i) satisfies the conditions 



H, 



2 3 



J 12 3 

then 



1 k 
1 i 



> 0, where k = 2, 3, 4, . . . and i = fe, fc+1, fc+2, . . . , (3.7) 



Pj-i(O)P^) = Pj-i( z ) 



ft 



2 3 
2 3 



Proof. Without loss of generality we assume j = 0. Suppose that 60 po = P-i- 

/2 3\ 
Then the minor /3q = i?o 00) nas two proportional rows, and is therefore zero. 

The converse we prove by contradiction. Let /3q = b\ — ai&o = and pa(z) ^ 
bop-i(z). Then there exists an integer i > 1 such that &j 7^ b$ai (since p-\{z) ^ 0). 
Therefore, according to (3.7) we have 



Consequently, 



H 



Ha 



2 3 4 

' 2 3 i + 2 



2 i + 1 



a aj 
60 fo 



>0. 



a ai a; 






6 &i 6^ 
a Oi_i 


= -ao 


a a, 
fro h 



<o, 



which contradicts the conditions (3.7). 



D 



Corollary 3.3. Let meromorphic functions p(z) and q(z) be such that p(0) = 1 
and g(0) Js 0. If the matrix H(p,q) satisfies the conditions (3.7), then F := 2 £ S. 
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Proof. For q(z) = this corollary is obvious. Suppose that q(z) ^ 0. Set 
po(z) '■= p(z) and p—i{z) '■— q{z) such that F(z) = Fq(z). 

Suppose that for some j ^ we have constructed the sequences p-i(z), . . . ,Pj{z) 
and Po, ■ ■ ■ , Pj-i > 0. By Lemma 2.1 the matrix Hj satisfies (3.7). Therefore, 
according to Lemma 3.2, we have two mutually exclusive possibilities: Pj > and 
Pj_i(0) • Pj(z) = pj-i(z). Consider the latter case. We have JT- ; = Pj{0) = 1 



and H 



(■!) 



H 



(4) 



• = with the notation (2.2). Additionally, the numbers Pi 

(k) 

are positive for all i = 0, 1, . . . , j — 1. By virtue of Lemma 2.1 the minors H a 
are positive for k = 2, . . . , j + 2, and zero for k > j + 2. So by Theorem 2.8 the 
considered function F(z) is a rational 5-function. 

If f3j > we can define the function pj + i(z) by (2.4). According to Lemma 2.1 
the matrix -ffj+i satisfies (3.7). So we can make the next step of this algorithm. 

(k) 

If this process is infinite, by Lemma 2.1 all the principal minors H are positive. 
Hence F eS. □ 



_ <M 

'-' """> '""" x ' v-i 
and q(z) = X)fc°=o bkZ k , normalized by the equality p(0) 
conditions are equivalent: 

(i) The infinite Hurwitz-type matrix 



Theorem 1.4. Consider the ratio F(z) — ^W of power series p(z) = X^^Lo a kZ k 

ao = 1. The following 



H(p,q) 



bo 


&1 


b 2 


b 3 


bi 


h ... 





«o 


a± 


(1-2 


a?, 


a^ 





bo 


h 


h 


b 3 


6 4 ... 








«() 


Oi 


«2 


a 3 ... 








bo 


b x 


b 2 


b 3 ... 



(3.8) 



is totally nonnegative. 
(ii) TTie matrix H{p, q) possesses the infinite factorization 

H(p,q)= lira (J(b ,Po) J(l,Px)--- J(l,Pj)) H(l,l)T(g), 



(3.9) 



converging in \\ ■ \\ p -norm for some p, < p ^ 1. .Here 60 S* and i/ie 
sequence {Pj)j^o is nonnegative, has a finite sum and contains no zeros 
followed by a nonzero entry, that is 

Po,Pi,...,P u -i>0, Pu = Pw+i = ■■■ =0, 



^ uj ^ 00, and 



J2 ft < °°- 
3=0 



The matrix T(g) denotes a totally nonnegative Toeplitz matrix of the 
form (1.3) with ones on its main diagonal. 
(iii) The ratio F(z) is a meromorphic S -function; its numerator q{z) and de- 
nominator p(z) are entire functions of genus up to a common meromor- 
phic factor g(z) of the form (1-4), g(0) = 1. 

Proof. We are proving as follows: (iii) ==>• (ii) => (i) => (iii). 

The implication (iii) =>■ (ii) is a consequence of Lemma 3.1 since total positivity 
of the matrix T(g) is provided by Theorem 1.3. 
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The factors in (3.9) are totally nonnegative and have at most a finite number of 
nonzero entries in each column. Therefore, the Cauchy-Binet formula is valid and 
the products of the form 

J(b , fa) J(l, /3i) • • • J(l, M H(l, 1) T(g), where < w < oo, 

are totally nonnegative. Moreover, the minors depend continuously on the matrix 
entries, so the entry-wise limit of totally nonnegative matrices is totally nonneg- 
ative itself. As a consequence, the matrix H (p, q) is totally nonnegative if the 
condition (ii) holds. That is the implication (ii) =>■ (i) is true. 

Let us show that (i) =>■ (iii). The matrices T(p) and T(q) are submatrices of 
H(p,q). Thus, if (i) is true, by Theorem 1.3 the series p(z) and q{z) converge to 
meromorphic functions of the form (1.4), so F(z) is a meromorphic function as well. 

Corollary 3.3 yields F £ S. So according to Lemma 2.3, the zeros (we denote 
their number by wi ^ oo) and poles (we denote their number by W2 ^ oo) of F(z) 
are real, simple and interlacing. Moreover, F(x) > F(0) ^ for x > 0, hence F(z) 
has only the zeros —t\, — T2, • • • , — t Ui and poles — u\, —u%, ■ ■ ■ > — ^u 2 , satisfying 
the following condition (cf. (2.18)) 

5^ n < ui < r 2 < <r 2 < t 3 < (3.10) 

Note that (3.10) implies the inequality coi — 1 ^ u>2 ^ 0J\. 

However, the functions p(z) and q(z) have the form (1.4), in particular they have 
no nonpositive poles. Therefore, all the numbers — n, — Ti, ■ ■ • , — t Wi are among 
the zeros of q(z), while —<j\, —02, ■••, — c W2 are among the zeros of p(z). As a 
result, 

q(z) = e' tlZ q{z) ) U- and p{z) = e^ z p{z) ) 

n,(i-t) ' n,(i-- 

where 71, 72, (av)„ and {Pp)^ are appropriately chosen positive numbers, 

ul1 ( z \ ^ 2 / z 
q(z) := b z j Y\_ ( H ) for j = 1 — sign 6 and p(z) := JJ H 

After cancellations in the fraction ——- we obtain Flz) — e lz , where 

p(z) p(z) 

7 : = 7i - 72- 

Let us show 7 = 0. Set 

G(z) := le-^F(z). 

For Wi,W2 < °° we can express the rational function G(z) as a sum of partial 
fractions and ascertain that it agrees with the expansion (2.6). So in this case G € S. 
Suppose that wi and 0J2 are infinite. With Arg : C \ {0} — > (— 7r, 7r] denoting the 
principal argument, we obtain the following for Imz > from (3.10) 

7T > Arg (to + z) > Arg (er + z) > Arg (n + z) > Arg (01 + z) > ■ ■ ■ > 0. 

Whenever ^ or to is positive, we obtain 

< Arg ( 1 + — ) - Arg ( 1 + — ) < Arg ( 1 + — ) - Arg ( 1 + -?— ) < vr. 
V W V °v) \ t v ) \ t v+ x) 



(3.11) 
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Therefore, if tq > 

OO -1 I Z OO -|i _2_ 

o < Arg [] y-±^ £ Ar § rri: = 

r/=0 en, i/=0 t^ 

|( Atg ( 1 + i)-A tg ( 1 + ±))< 

Arg | 1 + — ) - lim Arg { 1 + — ) < tt. (3.12) 

That is 

< Arg G(z) < tt when Imz > 0, (3.13) 



i.e. Ge5 since G(^) is real. If tq = we just replace all instances of ( 1 + — ) in 



TO 



inequalities (3.11) and (3.12) with z and obtain the same. This method to deduce 
the estimate (3.13) is taken from [5, Ch.IV §10, Lemma 11]. 

Now if 7 7^ 0, then Arg F(-ni/^) — Arg(—G(Tri/j)) and Ge5, which contradicts 
the inclusion F € S. Thus 7 = 0. □ 



4. Distribution of zeros and poles 
First, let us prove the following auxiliary fact. 

Lemma 4.1. Let p(z), q(z) be the formal power series p(z) = X^fco a k* and 
l( z ) = Sfc°=o bkZ k such that ao = and &o > 0. The matrix H(p,q) defined 
by (3.8) is totally nonnegative if and only if p{z) = and q(z) converges to a 
function of the form 



M = f ev ) ^, (4.1) 



where 7^0, a v , /3 M > for all n, v and J2„ ^7 + S M ~g~ < °° - 
Proof. If ao = then the total nonnegativity of H = -ff (p, g) implies 



< a, : = - — 
00 



ao a,; 
fro &i 



'M^+1, 1 '" W=LJS " 



so p(z) = 0. The Toeplitz matrix T(g) defined by (1.3) is totally nonnegative as a 
submatrix of H(p, q). Therefore, by Theorem 1.3, q{z) is of the form (4.1). 

Conversely, if p(z) = then any nonzero minor of H(p, q) is equal to a minor 
of T(q). According to Theorem 1.3 the matrix T(q) is totally nonnegative, hence 
H (p, q) is totally nonnegative as well. □ 

Theorem 1.2. A power series f(z) = X^fcLo A -Z ' c w ^ /o > converges to an 
entire function of the form 



1<U<U) 



M = foe» J] ( 1 + — )' ( 4 ' 2 
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2:; 



where 7 ^ 0, a v > for all v and Ylx<u<ui ^7" 
and only if the infinite matrix 



< 00 for some ui, $J u ^ 00, i/ 



2?/ 



//o /l 


h 


/a 


h ■ •■ 


/1 


2/2 


3/3 


4/4 ... 


/o 


h 


/2 


h ■■■ 





h 


2/2 


3/3 ... 



V : : 


/o 


h 


/a ■•■ 



is totally nonnegative. 

Proof. Let the matrix T>f = H(f',f) be totally nonnegative, where f'(z) denotes 
the formal derivative of f(z). If f\ = 0, by Lemma 4.1 f{z) = /o > 0, i.e. (4.2) is 
satisfied. 

Suppose that f\ ^ 0. Theorem 1.4 implies that /(z) and /'(z) converge in a 
neighbourhood of the origin. Moreover, for some meromorphic function g(z) of the 
form (4.1), the functions f(z) := ^rfj an d h := ^-4y are entire of genus 0, coprime 
and have only negative zeros. In particular, the poles of g(z) are positive (if any). 
Let us show g(z) has no poles. Observe that in the right-hand side of the expression 



h(z) = 



f(z)g'(z)+g(z)f'(z) 
9{z) 



f'(z) + f(z) 



9'{z) 
9{z) 



the logarithmic derivative of g(z) is multiplied by a function with no positive zeros. 
Therefore, each pole of g(z) must be a pole of h(z). But h(z) is an entire function, 
thus g(z) is entire and f(z) — f(z)g(z) can be represented as in (4.2). 

Conversely, let f{z) admit the representation (4.2). If f(z) is a constant then by 
Lemma 4.1 the matrix T>f is totally nonnegative. Suppose now that f{z) is not a 
constant and consider its logarithmic derivative 



F(z) 






7+ 



E — 



V 



^ w ^ 00. 



(4.3) 



x . The func- 

is an 



Each summand in the right-hand side of (4.3) is in 7?. , so F e 72." 
tion /(z) is non-constant, hence -F(z) ^ 0. Therefore, the functk 

"R.- function. Moreover, F{z) (and hence Wzt) has only negative poles and zeros. 
Consequently, j^p- is an 5-function. 

Since f{z) has the form (4.2), each common zero of f(z) and f'(z) is negative. 
In addition the functions e~ lz f(z) and e~ lz f'[z) are of genus 0. So the function 
fTzi = TTz) sa tisfies (iii) in Theorem 1.4. Consequently, the matrix 



■D f = H(f,f) = f 1 Hl^ 



is totally nonnegative. 



□ 
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H 



f 



V 



Let f(z) — J2T=o fk zk i /o > 0, be a real entire function. Define its infinite 
Hurwitz matrix 

I fa H fi fe fs 

o h h h h 

/o h h h 

o o h h h 

o O/0/2 h 



(4.4) 



For the minors of H. r we use the same notation as in Section 2, such that 



u 



H f 



2 3 
2 3 



Grommer in [8, §16, Satz IV] extended the Hurwitz criterion [11] to entire func- 
tions. However, he overlooked the condition on common zeros of odd and even parts 
(which was addressed by Krem in [14]). We only need the following particular case 
of this extension. 

Theorem 4.2 ([14, Theorem 12], [5, Ch.V §4]). Let a real entire function f(z), 
/(0) > 0, be of genus 1 or 0. Suppose that its even part (f(z) + /(— z))/2 and its 
odd part (f(z) — /(— z))/2 have no common zeros. 
Then the function f can be represented as 



m-oer- n (i+f) n (1 

1^/x^i V M/ 1<i/<w 2 V 



z 



1 



z 



(4.5) 



where $J Wi,u>2 ^ 00, 7^0 and C > 0, and its zeros satisfy the conditions 
x^ > for 1 ^ /j ^ wi, Rca„ > 0, Ima„ > for 1 ^ v ^ u>2, 



— < 00 and 



V Re — 



< CO, 



i/ and only if 



/(2)/ 



/(3)/ 



H^(/),H^(/),...,^ 1+2 " 2+1) (/)>0, 

t,(wi+2u;2+2) / j.s ^,(cji+2w 2 +3) / (-n, ... — (1 



Note that the restriction on the genus of f(z) implies the additional condition 



Ei 



^u^ui2 \a u 



< 00. Based on Theorems 1.4 and 4.2 we deduce the following fact. 



Theorem 1.1. Given a power series f(z) — z- 7 Efc°=o f kZ ^ > w ^ lere /o > and j is 
a nonnegative integer, the infinite matrix T-Lf defined by (4.4) is totally nonnegative 
if and only if the series f(z) converges to a function of the form 



f(z) = Cz 3 e' flZ+ ' f2Z 



iw< wi (1 + 1) ni««. ( l + t)( 1 + 1) 



for some ui\, u>2 and lo%, ^ u%, o>2; W3 ^ 00. Here C > 0, 

71,72 ^ 0, x v , y\ > 0, Rea^ ^ and Ima„ > for all /1, ^, A, 

. - -- > K<' 1 — i i- > - ' 

■I U 



f-i x u . Af 6 I a„ J . Af 



l^^^Wl 



l^y^Wa 



1<C f^W'2 



> — «■ < OO. 



(4.6) 

(4.7) 
(4.8) 



1^A^CJ 3 
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Proof. Suppose that f(z) is represented as (4.6). We can express it as 

f(z) = Cz J g(z 2 )h(z), where 

Z \ -n- ( ^ Z \ ( Z 



h(z)-.=e^ n i+ f n k i+ i- and ^ 






(- 2 )-^ n i+sfcr / n h-Sr • (4- 10 ) 






■Va 



Note that 5(2) has the form (4.1), so by Theorem 1.3 its Toeplitz matrix T(g) is 
totally nonnegative. 

Split h(z) into the odd part zh (z 2 ) and the even part h e (z 2 ) so that 

h(z) = h e (z 2 ) + zh (z 2 ). 

The function /i(z) is of genus not exceeding 1 as well as h e (z 2 ) and h (z 2 ). This im- 
plies that the genus of h (z) and h e (z) is 0. Indeed, for example, the function h e (z 2 ) 
with zeros ±<5„, n = 1,2,..., can be represented as the Weierstrafi product 



z 
- \ T ( * . z 



And since it depends only on z 2 , we necessarily have c = in this representation, 
which implies 

*.M-n(i-j| 

Let us show h (z) and h e (z) are coprime. Denote r := gcd(/i , /i e ) if h (z) ^ 0. If 
/i (z) = we set r(z) := h e (z). Assume that r(z) ^ 1. So it has zeros, since its genus 
is zero and r(0) = h e (0) — 1. However, if r(z ) = then h(^/zo) = h(—-y/zo) = 0. 
Since one of the points ±y/zo ^ is in the closed right half of the complex plane 
(independently of the branch of the square root) we get a contradiction to (4.9). 

If h (z) = then H(h ,h e ) = H(0,1) is totally nonnegative. This implies the 
total nonnegativity of the matrix 

U f = C H{gh ,gh e ) = C H(h ,h e )T(g). 

If h (z) e^ 0, then the function h(z) has the form (4.9) and its odd and even parts are 
coprime. That is, h(z) satisfies the conditions of Theorem 4.2. Therefore, % h = 
H( 2S) (h ,h e ), W h = H (3 )(h ,h e ),. . . is a positive sequence possibly followed by 
zeros. Thus, by Theorem 2.8, we obtain jf- £ S. Then the applying of Theorem 1.4 
to the function jf- gives the total nonnegativity of the matrices H{h ,h e ) and, 
consequently, Hj. 

Let us prove the converse. Suppose that the Hurwitz matrix T-Lf is totally non- 
negative. We can split the series f ~ 1 z~ : > f(z) = X^fclo f zk ln ^° the even part q(z 2 ) 
and the odd part zp(z 2 ) so that f(z) can be expressed as follows 

f(z) = f z j (q(z 2 ) + zp(z 2 )) . 

It gives Hf = f H(p,q). 
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If /i > 0, then since the matrix H(p,q) is totally nonnegative, by Theorem 1.4 
there exists a meromorphic function g(z) of the form (4.1) such that p{z) — ^l 

and q{z) = $¥l are coprime entire functions of genus 0. Moreover, the ratio ^l is 
an 5-function. Let 

f(z):=q(z 2 ) + zp(z 2 ). 

By Theorem 2.8 the minors %~ = H^ 2 \p,q), H~ — H^'{p,q),. . . form a pos- 
itive sequence possibly followed by zeros. Since q{z) and p[z) are coprime, the 
function f(z) has the form (4.5) by Theorem 4.2. 

If /i = then according to Lemma 4.1, p(z) = and q(z) has the form (4.1). 
Here we set g(z) := q(z) so that f(z) = 1. 

Now consider both cases /i = and f\ > 0. We showed that f(z) can be 
represented as in (4.5), while g(z) can be represented as in (4.1). That is, after 
the appropriate renaming of zeros and poles, the function f(z) has the form (4.9), 
the function g(z 2 ) has the form (4.10) and the conditions (4.7)-(4.8) are satisfied. 
Thereby 

f(z) = f z j g(z 2 )f(z) 
can be represented as in (4.6). □ 
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